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Abstract 

We analyze the radiative pion decay tt + — > e + is e j within nonlocal chiral quark models that 
include wave function renormalization. In this framework we calculate the vector and axial-vector 
form factors Fy and Fa at q 2 = — where q 2 is the e + v e squared invariant mass — and the slope a 
of Fv(q 2 ) at q 2 — > 0. The calculations are carried out considering different nonlocal form factors, 
in particular those taken from lattice QCD evaluations, showing a reasonable agreement with the 
corresponding experimental data. The comparison of our results with those obtained in the (local) 
NJL model and the relation of Fy and a with the form factor in 7r° — > 7*7 decays are discussed. 

PACS numbers: 12.39.Ki, 11.30.Rd, 13.20.Cz 
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The radiative pion decay 7r + — > e + v e r y is a very interesting process from different points 
of view. According to the standard description, the corresponding decay amplitude consists 
of the inner bremsstrahlung (IB) and structure-dependent (SD) terms. The IB contribution 
corresponds to the situation in which the photon is radiated by the electrically charged ex- 
ternal legs, either pion or lepton, while the SD terms are associated with the photon emission 
from intermediate states generated by strong interactions. The latter can be parameterized 
through the introduction of vector and axial- vector form factors, Fy{q 2 ) and Fa{q 2 ) respec- 
tively, where q 2 is the squared invariant mass of the e + z/ e pair Since 7r + — > e + u e is 
helicity suppressed, same happens to the IB contribution to its radiative counterpart, and, 
consequently, tt + — > e + v e r y turns out to be an appropriate channel to uncover the nonper- 
turbative SD amplitude. Recent measurements^ of the ir + — > e + v e r y branching ratio over a 
wide region of phase space yields F v (0) = 0.0258(17), F A {0) = 0.0117(17) and a = 0.10(6), 
where a is related to the dependence of F v on q 2 parameterized as 

F v {q 2 ) ~ F v (0) (1 + a q 2 /ml + ) for q 2 < m 2 w+ . (1) 

From the point of view of the physics of weak interactions, this determination has provided 
a further check of conserved-vector-current (CVC) hypothesis. In fact, the value of -FV(O) 

B'ven above is in good agreement to that extracted from the analysis of the 7r° — > 77 decay 
]. Moreover, as stressed in Ref. 2]], a good description of the data appears to be possible 
without the need to include extra tensor contributions that arise in several extensions of 
the Standard Model suggested in the literature (4, 5]. From the side of strong interaction 
physics, the pion radiative decays have been analyzed using Chiral Perturbation Theory jo] 
and effective meson lagrangian methods [7]. However, there still remains the question of 
how the associated form factors are connected to the underlying quark structure. Due to 
the nonperturbative nature of the quark-gluon interactions in the low-energy domain, to 
address this issue one is forced to deal with models that treat quark interactions in some 



effective way. Among these, the Nambu- 
an schematic effective theory for QCD 



ona-Lasinio (NJL) model has been widely used as 



lOj , allowing e.g. the description of light mesons 



as fermion-antifermion composite states. In the NJL model quarks interact through a local, 
chiral invariant four-fermion coupling. The corresponding predictions for the vector and 



axial- vector form factors at q 2 = have been calculated in Ref. ll|, yielding -FV(O) = 0.0244 



and Fa(0) = 0.0241. An extension of that calculation leads to a = 0.032. As we see, 
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while the predictions in the vector sector are in reasonable agreement with the measured 
values (the prediction for a being slightly below the empirical range), the calculated value 
of Fa{0) is a factor 2 too large. As a way to improve upon the NJL model, extensions which 
include nonlocal interactions have been proposed (see Ref. 12] and references therein). In 
fact, nonlocality arises naturally in several well established approaches to low energy quark 



dynamics. This is e.g. the case of the instanton liquid model [13 
resummation techniques [14], and also lattice QCD calculations 



and the Schwinger- Dyson 
indicate that quark 



interactions should act over a certain range in momentum space. Indeed, nonlocal chiral 



quark models have been successfully applied to study different hadron observables [18H26]. 
The aim of the present work is to investigate the predictions of this type of models for 
the measured quantities associated to vector and axial- vector form factors in tt + — > e + z/ e 7 
decays. 

In general, the amplitude for the process tt + — » e + v e (q) + 'y(k) can be written as [l| 
Gp 
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with (q + k) 2 = m 2 , k 2 = 0. Here Gp and 9c stand for the Fermi constant and the Cabibbo 



angle, respectively; is the photon polarization vector, Z M is the lepton current, L ail is 
a lepton tensor, and Fv(q 2 ) and FA{q 2 ) denote the vector and axial- vector hadronic form 
factors mentioned above. In this work we are interested in the study of the predictions 
for the measured quantities associated with these form factors in the context of nonlocal 
chiral models. We consider SU(2) chiral models that include wave function renormalization, 
defined by the following Euclidean action 24. 126] 



S T 



d 4 x <ijj(x)(—i^ + m c ) , ijj(x) 
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ja(x)j a (x) +jp(x)j P {x) 



(3) 



Here m c is the current quark mass, which is assumed to be equal for u and d quarks, while 
the nonlocal currents j a (x), jp(x) are given by 



Ja{X) 



J d 4 z g(z) ip (x+ r a ip (x - I 



3p( x ) = / <^ z f( z ) $ l x + 



^L^(x-- 

2x/l 2 



(4) 
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where T a = (l,ij 5 r) and u(x') d v(x) = u(x')d x v(x) — d x >u(x')v(x). The functions g(z) and 
f(z) in Eq. (jlj) are nonlocal covariant form factors characterizing the corresponding interac- 
tions. In what follows it is convenient to Fourier transform g(z) and f(z) into momentum 
space. Note that Lorentz invariance implies that the Fourier transforms g p and f p can only 
be functions of p 2 . 

In order to deal with meson degrees of freedom, one can perform a standard bosoniza- 
tion of the theory. This is done by considering the corresponding partition function 
Z = f VipUip exp[— Se], and introducing auxiliary fields a±(x), 02(2;), tt(x), where 01,2(2;) 
and tt(x) are scalar and pseudoscalar mesons, respectively. An effective bosonized action is 
obtained once the fermion fields are integrated out. To treat that bosonic action we assume, 
as customary, that a\^, fields have nontrivial translational invariant mean field values (7^2, 
while the mean field values of pseudoscalar fields 7Tj are zero. Thus we write 



a Ax) — u\ + 5a Ax) 



a 2 (x) = x p a 2 + 5a 2 (x) 



7r(x) = 5n(x) 



(5) 



Replacing in the bosonized effective action and expanding in powers of meson fluctuations 
we get 

cbos qMFA 1 riquad 



Here the mean field action per unit volume reads 



cMFA 

°E 

~ 2G 



1 (af + al) - AN C [ h, 



(27T) 



(6) 



with 



z p = (1- a 2 f p ) 1 , m p = z p (m c + a x g p ) . (7) 

The minimization of 5*^ FA with respect to a\, 2 leads to the corresponding gap equations. 
The quadratic terms can be written as 

1 f d 4 p 



^iquad 



J2 G M (p 2 ) 5M(p) 5M(-p) 



(8) 



M=0,0' ,7T 

where a and a' fields are meson mass eigenstates, defined in such a way that there is no 
a — a' mixing at the level of the quadratic action. The explicit expressions of Gm{p 2 ), as well 
as those of the gap equations mentioned above, can be found in Ref. 26J. Meson masses can 
be obtained by solving the equation Gm(— ht>m) — 0, while on-shell meson-quark coupling 
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constants gu q q are given by 

_ 2 dG M {p 2 



9Mqq 



dp 

As in Ref. 



(9) 



26], we will consider here different functional dependencies for the form factors 
g p and f p . First, we consider a relatively simple case in which there is no wave function 
renormalization of the quark propagator, i.e. f p = 0, z p = 1, and we take an often used 
exponential parameterization for g p , 

g p = ex V (-p 2 /A 2 ) . (10) 

The model parameters m c , Gs and A are determined by fitting the pion mass and decay 
constant to their empirical values = 139 MeV and f n = 92.4 MeV, and fixing the chiral 
condensate to the phenomenologically acceptable value {qq) 1 ' 3 = —240 MeV. In what follows 
we refer to this choice of model parameters as Set A. Second, we consider a more general 
case that includes the wave function renormalization of the quark propagator. We keep the 
exponential shape ([10]) for the form factor g p and assume also an exponential form for f p , 
namely 

j p = exp (-pVA?) . (11) 

Note that the range (in momentum space) of the nonlocality in each channel is determined 
by the parameters A and A 1; respectively. As in the previous case, model parameters are 
determined so as to reproduce the desired values of m^, f n and (qq) 1 ^ 3 - The form factor 
f p introduces now an additional free parameter A 1; consequently we introduce as a fourth 
requirement the condition z p (0) = 0.7, which is within the range of values suggested by 
recent lattice calculations [8, 10]. This choice of model parameters and form factors will be 
referred to as parameterization Set B. Finally, we consider a different functional form for 
the form factors, given by 

_ 1 + a z a m fmip) ~ m a z f z {p) j. _ 1 + a z j. ,s ^ 

9p l + a z f z (p) a m -ma z p 1 + a z f z {p) zP 



where 



Up) = i + ( P 2 /Kf 2 • m = [i + (p 2 /m~ b, ' z ■ (13) 
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As shown in Ref. 



26|, taking m c = 2.37 MeV, a m = 309 MeV, a z = -0.3, A = 850 MeV 



and Ai = 1400 MeV one can very well reproduce the momentum dependence of mass and 



renormalization functions obtained in lattice calculations, as well as the physical values of 
and f n . In what follows we will refer to this choice of model parameters as parameterization 
Set C. The parameter values for all three parameter sets, as well as the corresponding 



26]. 



predictions for several meson properties, can be found in Ref. 

In order to derive the form factors we are interested in, one should "gauge" the effective 
action Se by introducing the electromagnetic field A^(x) and the charged weak fields W^(x). 
For a local theory this "gauging" procedure is usually done by performing the replacement 

9^9 M +iG„(x), (14) 

where 

W) = 2 I 3 JM x )+9w— ^ , (15) 

with <jrjy/(8Mjp) = Gp cos 6 C / \/2 and r 1 * 1 = (r 1 ± ir 2 )/2. In the present case — owing to 
the nonlocality of the involved fields — one has to perform additional replacements in the 
interaction terms, namely 

i/>(x-z/2) ->■ W G (x,x- z/2) ip{x-z/2), 

^{x + z/2) -> ^(x + z/2) W G (x + z/2,x) . (16) 

Here x and z are the variables appearing in the definitions of the nonlocal currents [see 
Eq.fll])], and the function W G (x,y) is defined by 



W G (x,y) = P exp 



y 

X 



(17) 



where r runs over an arbitrary path connecting x with y. 

Once the gauged effective action is built, the explicit expressions for the vector and 
axial-vector form factors can be obtained by expanding to leading order in the product 
5tt + A^W„ . This is a rather lengthy calculation that can be simplified by considering 
q 2 rn 2 + , as needed to make predictions for the measured quantities F v (0), Fa(0) and a. 

The vector form factor is obtained from the triangle diagram represented in Fig. la. As 
stated, Fv(q 2 ) can be expanded at leading order in q 2 as in Eq. (JTJ, just changing q 2 — > —q 2 
since we are working in Euclidean space. Performing such an expansion we obtain 
F v (0) y/2g wq qN c f dH (z e + z iic )(z e + z tq ) \ £ 2 (m lQ - m e m lfC -m t 



(2tt) 4 9i0 D t D iK D lQ Zi 



2 V. K 



18) 



where K and Q are the photon and e + u e pair momenta, respectively, taken at q 2 = 0. 
In Eq. ffT8l) we have used the definitions Di = £ 2 + mj, £q = £ + Q, £k — £ — K and 
A) — (^q + ^x)/2. For convenience we also define £q = £ ± Q/2, £^ = £ ± K/2 and 
£± Q = £±(K + Q) /2, which will be used below. 

The calculation of a is somewhat more involved. In particular, the result depends on the 
integration path appearing in the nonlocal contribution to the quark— gauge boson vertices 
[see Eq. ffTTl) ]. One obtains: 



d 4 £ (z t + z tK ) {z t + z e ) 

9e — FT~R ^ B(£,K,Q) 



(27T) 
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(20) 



Here primes stand for derivatives, e.g. g\ = dgi/d£ 2 . The functions ct^Q and fti t Q are, in 
general, path dependent. Here, for simplicity, we choose to use the "straight line path" for 
which they read 



®i,q = 0i J d\ A g'i+ + x q , 0e,Q = &2 J d\ A f' e + + x 



(21) 



The axial-vector form factor receives not only a contribution from the triangle diagram 
in Fig. la (as occurs in the local NJL model) but also from other diagrams, which are 
represented in Figs.lb-le. Thus, for q 2 = one has 



Fa(0) = j> A (0)| t 



(22) 



The contribution from the triangle diagram is given by 



^a(O) 
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where 
Aa =2(4 



K -i Q -i 
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Q 
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ze 



while the remaining contributions are given by 
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Fa(0) 
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where 
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, K-(e + lQ) 

Note that, contrary to what happens with Fy(0), the axial-vector form factor depends on 
the path even at q 2 = 0. This is due to the contribution of the diagram of Fig. le. The 
expression for j(£, K, Q) given above corresponds to the "straight line path" choice. 

Before discussing our predictions for the form factors associated with the charged pion 
radiative weak decay, let us note that one can also consider the related decay processes 
7T° — > 77 and 7r° — > e + e~7. The amplitude of these processes contains the 7T°77 vertex form 
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factor F 71 " 77 * (g 2 ), where q 2 is now the invariant mass of the virtual photon. As in the case 
of Fy(q 2 ) in the decay 7r + — > e+w-f, one can perform an expansion for low q 2 : 

p^r „ ^77* ( ) (! + a > g 2 / m 2 ) ; W ith q 2 < m 2 . (28) 

Experimental measurements lead to i^ 77 *(0) = 0.284(8) GeV^ 1 and a' = 0.032(4) fl. 
Here it is important to recall that, in the chiral limit, the anomaly leads to the constraint 
F' r77 *(0) = l/(47r 2 / 7r ) ~ 0.274 GeV -1 . Moreover, using CVC and isospin arguments it is 
easy to prove that F v (0) = F 71 ' 77 * (0) / y2. We stress that our model satisfies this relation 
as well as the anomaly constraint. In addition, it is easy to see that it leads to an analytical 
expression for a' which coincides with that given in Eq. f fl9|) for a. 

We discuss now our numerical results. In Table I we list the predicted values of F v (0) 
and the slope a for the three parameterizations considered in this work. We also include the 
available empirical data, as well as the (local) NJL model predictions. In the case of Fy(0), 
we observe that the predictions of all three parameterizations of the nonlocal model are in 
good agreement with the empirical value and with the value obtained in the NJL model. 
This is hardly surprising, given the chiral limit constraints mentioned above. Regarding the 
slope parameter we find a too small value for all three parameterizations. There is a small 
improvement by going from Set A to Set B (which implies the introduction of the wave 
function renormalization of the quark propagator), and a larger improvement is obtained 
by going from set B to C. One of the peculiarities of the nonlocal model is that the results 
can be dependent on the path used in Eq. (flTl) . As usual, we have chosen a straight line 
path for the calculations, obtaining the expressions given in Eq. fl2TT) for the path dependent 
quantities ae t Q and /^q. To gauge the importance of this path dependence we have evaluated 
the contribution of the corresponding terms, obtaining that it represents less than 3% for 
parameterizations A and B and less than 0.5% for Set C. As for the predictions for a' we 
note that, even though in our model its analytical expression coincides with that for a, the 
corresponding numerical values are somewhat different since in the case of a' we evaluate 
Eq. (HHD at P 2 — tti 2 o , whereas for a it is evaluated at P 2 — m 2 + . This dependence on 
the mass of the pion reduces the value of a' in comparison with that of a by about 10% for 
Sets A and B and 5% for Set C. We observe in Table I that a' is better reproduced than 
a, but the large errors in the experimental determination of the latter prevent us to take 
definite conclusions. In fact, the discrepancies between our predictions for a and a' and their 
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experimental values are not unexpected, since no vector interaction has been included in our 



model. For example, using the results of Ref. |28j we obtain that the vector contribution 
to the 7T° — > 77* process has the right order of magnitude and sign needed to account for 
these discrepancies. Here it is interesting to point out that the contribution from the vector 
channel could be different for a and a', even preserving the isospin symmetry. This would be 
achieved if one has different interactions in the vector-isoscalar channel and in the vector- 
isovector channel. On the other hand, it can be seen that the NJL model is able to provide 
a very good prediction for a'. However, it is worth to notice that in the NJL framework this 
value is essentially given by the relation a' ~ {12m 2 q ) [27J, and the value of m q (dressed 
quark mass) turns out to be quite dependent on model parameterizations. The NJL values 
quoted here correspond to a Pauli-Villars regularization. 

In Table II we give our numerical results for Fa(0), quoting the contribution of each 
diagram. We have also quoted in Table II the result obtained in the NJL model, where the 
triangle diagram is the only one that survives. We observe that although several diagrams 
contribute in the case of the nonlocal models, the triangle diagram turns out to be the 
dominant one. The contribution coming from diagrams (b) to (e) is less than 2% of the total 
result, in all three parameterizations. In particular, diagram (e), the only one with a path 
dependent term, contributes with less than 1% of the full result. It is interesting to mention 
that Fa(0) has also been evaluated in the spectral quark model, obtaining Fa(0) ~ F v (0) 



as in the NJL model 29|. In this sense, the fact that nonlocal models of the type considered 
here lead to values of Fa{0) which are significantly different from those of -FV(O) appears as 
a quite important result. A similar conclusion has been obtained in a simplified calculation 
in the chiral limit 30] . Given the triangle diagram dominance mentioned above, the origin 



of this difference can be traced back to the different dressing of the 7 M and 7^75 terms in the 



coupling of the W to the quarks [26|. Comparing our prediction for Fa(0) with the empirical 
value, we observe that the introduction of the nonlocality gives half of the difference needed 
in the case of Set C and is near to exhaust this difference for Set A, which gives the best 
result. In any case, as in the case of the slope of the vector form factor, additional effects 
that have been neglected in the present calculation might help in improving the agreement 
with the empirical value. 
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TABLE I: Results for Fy(0) and a. All results should be multiplied by 10~ 2 . In column 5 we give 
the empirical values of F^ 11 * (0)/\/2 and a'. Note that in our model a = a'. In column 6 we 
give the prediction of the local Nambu-Jona-Lasinio (NJL) model. 





Set A 


Set B 


Set C 


Exp [2] 


Exp (vr° -)■ 77*) [3] 


NJL 


F v (0) 


2.697 


2.693 


2.695 


2.58 ± 0.17 


2.80±0.08 


2.441 


a 


1.651 


1.726 


2.011 


10 ± 6 


3.2±0.4 


3.244 



TABLE II: Results for F A (0). All results should be multiplied by 10 2 . In column 5 we give the 
prediction of the local Nambu-Jona-Lasinio (NJL) model. 





Set A 


Set B 


Set C 


Exp [2] 


NJL 


F A (0)\a 


1.300 


1.591 


1.804 




2.409 


F A (0)\b 


0.067 


0.047 


0.031 






F A (0)\c 


-0.0002 


-0.0002 


-0.0001 






F A (0)\ d 


-0.044 


-0.036 


-0.026 






F A (0)\e 


-0.003 


0.013 


0.017 






F A (0) 


1.319 


1.614 


1.825 


1.19 ± 0.01 


2.409 
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(a) 




FIG. 1: Diagrammatic representation of the possible contributions to ir + — > e + v e ^ decay. Double- 
dashed lines, wavy lines and single-dashed lines represent the e v e pair, the outgoing photon and 
the decaying pion, respectively. While for the vector form factor only the contribution from the 
triangle diagram (a) is nonvanishing, in the case of the axial-vector form factor all five diagrams 
contribute. 
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